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ABSTRACT 

We present a study of the angular correlation function as measured in the William 
Herschel Deep Field, a high galactic latitude field which has been the subject of an ex- 
tensive observing campaign from optical to infrared wavelengths. It covers 50 arcmin^ 
and with it we are able to investigate the scaling of the angular correlation function 
to i? ~ 28, i?, J ~ 26, K ~ 20 and H ~ 22.5. We compare our measurements to re- 
sults obtained from the smaller Hubble Deep Field. To interpret our results, we use a 
model which correctly predicts colours, number counts and redshift distributions for 
the faint galaxy population. We find that at fixed separation the amplitude of o-!{9) 
measured in BRI bandpasses is lower than the predictions of a model containing with 
no luminosity evolution and stable clustering growth in proper co-ordinates. However, 
in the near-infrared bandpasses, our measurements are consistent with the predictions 
of an essentially non-evolving K— selected galaxy redshift distribution. In the range 
i? ~ 27 — 28 we find that our correlation amplitudes are independent of magnitude, 
which is consistent with the observed flattening of the number count slope and cor- 
respondingly slower increase of the cosmological volume element expected at these 
magnitudes. 

If our luminosity evolution models provide a correct description of the underlying 
redshift distributions (and comparisons to available observations at brighter magni- 
tudes suggest they do), then our measurements in all bandpasses are consistent with 
a rapid growth of galaxy clustering (0 < e < 2 in the normal parametrisation) on 
the sub-Mpc scales which our survey probes. We demonstrate that this rapid growth 
of clustering is consistent with the predictions of biased models of galaxy formation, 
which indicate that a rapid rate of clustering growth is expected for the intrinsically 
faint galaxies which dominate our survey. 



1 INTRODUCTION 



The projected two-point galaxy correlation function io{9) 
has proved to be one of the most enduring statistics in ob- 
servational cosmology. This is a consequence of the relative 
ease with which it may be measured; for each galaxy, all 
one requires is positions and magnitudes. Starting with the 
early st udies of clustering in the local universe using Schmidt 
plates ( Groth fc Peebles 1977) to more recent works using 
CCD-based detectors (Efstathiou ct al. 1991; t^ritchet fc In- 



cosmology and knowledge of the underlying redshift distri- 
butions for each magnitude-limited sample. 

In this paper we will investigate the projected angu- 
lar clustering of the faint field galaxy population. We char- 
acterise galaxy clustering as a function of sample limiting 
magnitude in BRIKH bandpasses. Our primary dataset is 
a deep, ground based survey of an area called 'the William 
Herschel deep field' ( WHDF). This has been described in 



fante 1 )92) these studies have probed galaxy clustering to 
very faint magnitudes. Normally, these surveys measure how 
the amplitude of the projected angular correlation function 
at a fixed angular separation, A^^, varies as a function of sam- 
ple limiting magnitude - the "scaling relation" . Usually, this 
relation has been parametrised in terms of "epsilon mod- 
els" in which the three-dimensional correlation length rpfz 



sever al recent papers (Metcalfe et al. 1996; McCracken et al 
200C ) . Covering ~ 50 arcmin'^ this survey comprises an area 
~ 10 times larger area than the separate HDF-N and HDF- 
S fields. For comparison, we also present a complementary 
analysis of clustering amplitudes measured in these smaller 
fields, utilising the catalogues produced in Metcalfe et al 
(2000). Although similar studi e s of uj{9) exist in the lit- 



m 



scales I lonotonically with redshift (ICroth fc Peebles I977I erd, Small, fc Mould 19941 ; |Hudon fc Lilly 1996 



erature (Efstathiou et al. 1991 



Phillipps et al. 197q). These models also require a choice of 



Roche et al. 1993; prain 



Woods 



Fahlman 1997) our survey differs primarily in its depth 
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{B ~ 28) and broad wavelength coverage (in this analysis 
we consider samples selected in BRIK bandpasses). 

To interpret our results we use redshift distributions 
derived from the luminosity ev olution models we have de- 
scribed in our pre vious papers ((McCracken et al. 2000; Met- 
calfe et al. 1996). These models are able to reproduce all 



the observable quantities of the faint field galaxy popula- 
tion (counts, colours, redshift distributions), at least for low 
Qo universes and within current observational uncertainties 



(Metcalfe et al. 1996); it is these successes which give us con- 



fidence in using our models as probe of the clustering his- 
tory of the Universe, rather than using our measurements 
of Aui as a probe of the underlying redshift distributions. 
In our models, high Slo Universes can be accommodated by 
the model if we add an extra population of low luminosity 
galaxies with constant star-formation rates which boost th e 
counts at faint {B > 25™) magnitude levels (Campos 1997). 
We als o consider flat cosmologies with A ^ 0. For reference. 



the sea ing relation computed for a model with stable clus- 
tering and no luminosity evolution is also presented. 

Models such as those presented in this paper are rela- 
tively successful in describing clustering measurements per- 
formed on deep blank-field surveys like the one detailed in 
this work ( [Roche et al. 1993 ; Brainerd, Small, & Mould 



1994|)~]h owever, observat ions of the cluster ing properties 



of Lyman-break galaxies (Madau et al. 1996) indicate that 
these objects hav e comparable clustering properties (Gi- 
avaliscc et al. 1998) to some classes of locally observed galax- 
ies, making such objects initially difficult to understand in 
terms of this monotonic scaling of ro with redshift. We will 
explain how these observations can be understood in the 
context of the results presented in this paper. 

Our paper is organised as follows: in Section ^ we de- 
scribe in outline the preparation of our datasets; in Section ^ 
we describe the techniques we use to measure and analyse 
our data; in Section |^ we present our measurements of the 
projected correlation function in flve bandpasses in compar- 
ison with previous work and investigate if our errors esti- 
mates are realistic; in Section ^ we compare our correlation 
measurements with the predictions of our evolutionary mod- 
els; and finally, in Section ^ we outline the main conclusions 
from this work. 



2 OBSERVATIONS AND CATALOGUES 

Full details of the optical observations comprising the 
WHDF will be presented in a forthcoming paper (Metcalfe 
et al 2000). A subset of our infrared observa tions of the 
WHDF is described in 



McCracken et al. (2000) which com 



prises the K < 20 UKIRT observations. Additional infrared 
observations at Calar Alto Observatory produced a second 
catalogue limited at H < 22.5 which will be fully described 
in a separate paper (McCracken et al, in preparation). In 
this section we will briefly describe our object detection and 
photometry techniques which are very similar to that used 



in our previous galaxy counts papers (Metcalfe et al. 1991 



Jones et al. 1991; Metcalfe et al. 1995; McCracken et al 



200C ). All our optical data discussed in this paper was taken 
at the William Herschel Telescope (WHT), with the excep- 
tion of a short /—band exposure made at the Isaac Newton 
Telescope (INT). 



After bias subtraction and flat-fielding, the sky back- 
ground is removed and isophotal image detection is carried 
out. These images are then removed from the frame, replaced 
by a local sky value, and the resulting frame smoothed heav- 
ily before being subtracted from the original. This produces 
a very flat backg round. The isophotal detection process is 
then repeated. A Kron (198C)-type pseudo-total magnitude 



is then calculated for each image, using a local value of sky. 

Table ^ shows the magnitude limits for our fields. As in 
our previous papers the minimum Kron radius is set to be 
that for an unresolved image of high signal-to-noise, and the 
correction to total is the light outside this minimum radius 
for such an image. Our measurement limits give the total 
magnitudes of unresolved objects which are a Scr detection 
inside the minimum radius (which is typically ~ 1.4" for 
the WHDF data). Star-galaxy separation was done on the 
B frame using the difference between the t otal magnitude 
and th at inside a 1" aperture, as described in Metcalfe et al. 
(1991). This enabled us to separate to B ^ 24™. Some addi- 
tional very red stars were identified from the R and I frames. 
As the WHDF is at high galactic latitude the stellar con- 
tamination should in any case be quite low. For the purposes 
of measuring the correlation function, masks files were also 
constructed to cover regions containing bright galaxies or 
stars. The area of the field affected by such bright objects is 
less than 10% of the total. 

Similar methods were also used to generate catalogues 
from the north and south Hubble deep fields (i.e., we do 
not use any of the existing HDF catalogues but use our 
own independently written object detection and photometry 
software). One significant difference between the HDF data 
and our ground-based data is of course their much higher 
resolution. As described fully in Metcalfe et al (2000, in 
preparation), we visually inspect all detections on our HDF 
N/S data in an attempt to reduce the number of spurious 
entries in our catalogues. We also carry out a 'reassembly' 
process in which multiple detections on an individual galaxy 
are combined to produce a single detection. This admittedly 
subjective procedure is unavoidable in the HDF catalogues 
given the extremely high resolution of the data. 



3 METHODS AND TECHNIQUES 

3.1 Determining the angular correlation function 



We use the normal estimator of Landy & Szalay (1993), 
given in equation |l|. Here we follow the usual notation 
where DD indicates the number of galaxy-galaxy pairs, DR 
galaxy-random pairs and RR random-random pairs for a 
given angular separation and bin width: 



;(e) = 



DD - 2DR + RR 
RR 



(1) 



We find that, for a given survey sample, amplitudes 
measured by this e stimator are ve ry similar to those 
computed with the Hamilton (1993) estimator, io{9) = 
(DDRR/DR^) - 1. The DD/DR - 1 estimator, as used by 
Roche et al. (199^ ), gives consistently higher (at the 20—30% 



level) values for uj{9), over all our bins. This has been found 
by other authors and is i ndicative of the known biases in- 
heren t in this estimator ( Hamilton 1993| ; Landy & Szalay 
199^). 
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Filter U B R I K H 

Limit (So-) 26.8 27.9 26.3 23.5/25.6 20.0 22.5 

Area (arcmin^) 48.5 48.5 48.5 88.3/53.0 47.2 50 

Table 1. Photometric limits of the WHDF. The two magnitude limits in I refer to two separate surveys, one carried out at the INT 
(and covering 88.3 arcmin^ and the other based on WHT data. 



For a range of magnitude-limited samples of each cata- 
logue, Lo{0) is computed using equation |l| for a series of bins 
spaced in increments of 0.2 in log{6), where 9 is in degrees. 
As we have only observed one field we cannot use the field- 
to- field variance to estimate the errors in each bin; instead 
we implement a bootstra p-resampling technique ([Barrow, 
Sonod 4 fc Bhavsar 1984| ; [Ling, Barrow, fc Frenk 19861 ). In 
this method, the error in each bin is computed from the 
variance of the estimator as applied to a large (~ 200) num- 
ber of bootstrap-resampled catalogues. As expected, these 
bootstrap errors are larger (normally ~ x 2) than the normal 
\/7V Poisson counting errors. 

To allow comparison with other workers, we fit our mea- 
sured correlations as a function of angular separation to an 
expression of the form 



iu(e) = A^{e'^ - c) 



(2) 



where A^i is the amplitude of io{9) at 1° and C is the "in- 
tegral constraint" term. This term is a bias which arises 
because we are using each catalogue to determine the mean 
galaxy density and is particularly significant in our work be- 
cause the area which we survey is quite small. To calculate 
the integral constraint we use the expression 



C 



1 

^2 



(3) 



where 9 is the angular separation of each galaxy pair and 
dSli and dD.2 the solid angle subtended by each pair. If we 
assume a power-law correlation function, ui(d) oc 0-°-* we 
may calculate this quantity for our fields by direct integra- 
tion. Typically we find C ~ 13 for the WHDF and ~ 40 for 
the HDF (we must assume a slope for power-law correlation 
function as we cannot calculate it directly from this data; 
—0.8 allows us to compare our work with similar studies in 
the literature). 

The error on A^, the overall fit, is determined from the 



method of Marquardt (1963), as described in Press et al 
(19860 This method combines errors on each bin in an in- 
dependent manner to calculate the total error of the fit. 
Figure |l| shows fits made for the B— band catalogue. 

We determine correlation amplitudes for the Hubble 
deep field data using a similar procedure. In this case we 
fit our final power law to an average of the correlation func- 
tion determined independently on each of the three WFPC2 
chips. For our NICMOS correlation amplitude, we compute 
our correlation functions from the total numbers of pairs 
from both surveys. For all these space-based data sets the 
field of view is extremely small, and consequently the re- 
quired integral constraint correction is very large. Addition- 
ally, the small numbers of pairs involved means that fits are 
generally dependent on three or fewer bins, and for this rea- 
son our resulting correlation amplitudes determined from 
these data should be regarded as upper limits on the fitted 
amplitudes, rather than definitive measurements. In order to 
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Figure 1. u)(9) as measured for samples limited at iJ < 27"* and 
B < 28™ . The solid line shows the fit to ui{e) = A„(e-°'® - C) 
where C is the "integral constraint" term described in the text 
and A^j is the value of ij^{8) at 1° 



try to reduce problems from "merged" objects as described 
in Section ti we carry out our fits at angular separation > 1". 



3.2 Modelling the correlation function 

We would like to compare our measured correlation ampli- 
tudes with those of model predictions. In order to do this 
we must assume a functional form for the s patial correlation 
function. From the results of large surveys ( Groth fc Peebles 



1977 



Davis fc Peebles 1983; Maddox et al. 19901 ) it is found 



that ^(r) (the spatial correlation function) is well approxi- 
mated by ^(r) = {ro/r)"', at least for scales < 20h~^ Mpc. 
Projecting a model for ^(r) onto the two-dimensional distri- 
bution of galaxies measured by 1^(9) involves integrat ing this 
function over redshift space using Limber's formula ( Limber 



1953) 



We must parametrise the scalin g of the correlation func- 



tion 



with redshift. Early papers ( [Groth fc Peebles 197 



Phillipps ot al. 197J ) assumed a scaling of the form 



C(r,z)=M^)(^)' 



(4) 
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where 
h{z) = + 



(5) 



(in this case r is the proper distance); this relation has been 
used in many recent observationaUy-motivated studies inves- 



1991; |l 


,oche et al. 19931; 


[nfante & Pritchet 1995a 



Braincrd, Small, fc Mould 1994t 
Small 1998| 



Brainerd 



To derive an expression for lu{6), the projected corre- 
lation function, we note that for small angles, the rela tion 
between uj{9) and ^(r) becomes (Efstathiou et al. 1991) 

-r[(7-l)/2] A , 



617- 



(6) 



r(7/2) 

where F is the incomplete gamma function, 6 is the angular 
separation and A is given by 



A = 



where 



dz/ 



hjz) 

dY^{z) {dr{z)/dz 



(7) 



(8) 



where dA{z) is the angular diameter distance and dr{z)/dz 
is the derivative of the proper distance. 

Analysis of the aforementioned large local redshift sur- 
veys suggests that 7 — 1.8, leading to three cases of interest 
to us: clustering fixed in proper coordinates, in which case 
e = 0.0; clustering fixed in co-moving coordinates which 
gives e = —1.2. Finally, the predictions of linear theory give 
e — 1.0. This formalism has been widely used in many papers 
which invest igate the clustering of fain t field galaxies: see, 
for example. Infante & Pritchet (1995b), [Woods fc Fahlman 
(1997|] 

As we have already noted, in these "epsilon models" 
characterised by equation ^ the co-moving galaxy correla- 
tion length decreases monotonically with redshift (providing 
of course e > —1.2, which produces models with clustering 
constant in co-moving co-ordinates) However, several recent 
works lliave indicated that ttus may not be a realistic as- 



sumption" la th(^6t-(^titia.l athdifes, both N-b6dy tjihlulAttohs 



( IColin et al. 1999|) and semi- analytic models ([Baugh et al 



1999 |; [Ffauffmann et al. 1999| ) indicate that the co-movmg 
galaxy correlation length decreases until 2 ~ 1 — 2 after 



which i t increas es again. These theoretical studies (Gover- 
nato et al. 1998 ) also allow us to explain the high cluste ring 
amplit udes observed for Lyman break gala xies at z ~ 3 (Gi- 



avaliscc et al. 1998; Adelberger et al. 1998) as a consequence 



of their formation in highly biased environments. Further- 
more, the clustering growth is expected to be more rapid 
for less massive objects and and for clusterin g amplitudes 
measured on smaller scales (Baugh et al. 1999). 



Motivated by these works we also model our correla- 
tion amplitudes using a modification of equation ^. In place 
of the normal epsilon parametrisation, we have used in the 
relativistic Limber's equation a more general form for the 
evolution of ^(r, z), namely 



il + z) 



(9) 



where r-o°™(z) is the comoving correlation length at z. Thus, 
we have used 
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Figure 2. Normalised co-moving correlation length rQ°'^{z) as 
fitted to haloes of circular velocity V > 120 km~^ as identified 
by Colin et al in the large N-body simulation of Kravtsov and 
Klypin (1999). This line is closest to the predictions of e = 
model described in the text. 



hiz) 



\z) 



{l + z)ro 



(10) 



To illustrate the possible effect of modelling more exactly 
the evolution of the correlation function, we have used the 



evolution seen i n the large N-body simulation of Kravtsov 



& Klypin (1999); the semi-analytic models mentioned above 
produce a similar form for the evolution of ^(r, z) in their 
simulations. As our field sample is dominated by spirals, 
we have therefore considered the haloes of the simulation 
having velocity v > 120km~^. Also, as uj{9) for these deep 
fields has, as usual, been fitted to a —0.8 power law, we 
have converted the Colin et al. data to provide the same 
correlation strength as a —1.8 power law for ^{r,z) at a 
comoving separation of 0.3 Mpc, which at the depths of 
our data here corresponds roughly to the angular scale of our 
estimates for uj{9). Finally, to obtain the function, ro°'^{z), 
a spline fit was made to the converted Colin et al. data 
points with a simple linear extrapolation to redshifts larger 
than the maximum redshift, z = 5, for which they have 
estimated the correlation function for their simulation. Fig. 2 
plots the resulting form of the evolution used for ro°™{z) 
normalised to ro. In using this in Limber's equation, we have 
taken, as with Roche et al. (1993), ro — 4.3 Mpc, which 
is little different from the converted Colin et al. value of 
4.2 h-^ Mpc. 
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3.3 Calculating dn/dz 

From equation ^ we see that the amphtude of lo{9) depends 
on the redshift distribution, dn/dz. To produce these red- 
shift distributions we employ a pure luminosity evolution 
(PLE) model in which star-formation increases exponen- 
tially with look-back time. Earlier versions of these models 



are disnussed in um pieviuus papers (MeLcalfe el al. 1991 



Metcalje-e t al. lOOS j ; [Metcalfe et al. 199^ ), and a full deacrip 



with measurements from the literature (open symbols) . The 
solid line shows the predictions of the stable clustering, e = 
non-evolving (i.e., no luminosity evolution) model , com- 
puted assuming ro = 4.3/i~^ Mpc and qo — 0.05 (This value 
of ro was chosen to produce the correct clustering ampli- 
tude at brighte r magnitudes as measured from early Schmidt 
plate surveys (Jones, Shanks, & Fong 1987 



Stevenson et al 



1985). We adopt the same value of ro for all bandpasses; in 



tion of the model used in this paper is given in (McCracken 
et al. 2( 00). In this paper we assume Hq — 50 kms~ Mpc~ , 



Section 5.4 we discuss if this is an appropriate assumption 



although changing the value of Ho does not markedly af- 
fect our conclusions. Two values of the deceleration param- 
eter qo = 0.05 and qo = 0.5, are adopted, corresponding to 
open and fiat cosmologies respectively. The input parame- 
ters to our models consist of observed local galaxy param- 
eters (namely, rest-frame colours and luminosity functions) 
for each of the five morphological types (E/SO, Sab, Sbc, 
Scd and Sdm) we consider in our models. These morpho- 
logical types are divided into ear1y-typp (F,/SO/Sah) and 



spiral ( Bbc/Scd/Sdml and these two classes are each given 
a separate star-formation history, parametrised in terms of 
an e-folding time r. We compute the k + e corrections us- 
ing the models of Bruzual fc Chariot (1991^ ). We could, in 
principle, sub-divide the spirals into different morphologi- 
cal types each with different star formation histories but for 
simplicity we do not; [k 4- e) corrections for the different 
types are fairly similar to each other in these models in any 
case. Instead, taking a Sbc model as representative of all 



for our data.) 

Starting with the B— band, we note that here our 
WHDF sample reaches extremely high galaxy surface 
density — approaching ~ 10® gal deg~^ at _B = 28™, and 
furthermore it probes to the highest redshift; our low-go 
evolutionary models indicate that by B ~ 28 we reaching 
2med ~ 2. Morcovcr, our measurements of the B— band cor- 
relation function are significantly deeper than any previously 
published work. Our brightest bin, at B < 27.0, is in agree- 
ment with t he correlation amplitude measured by Metcalfe 
et al. (1995 ). Faintwards of B = 27, our correlation am- 



plitudes remain flat. The errors on our fltted correlations 
in B are relatively low in comparison with our other band- 
passes because at B < 28 we detect ~ 6000 galaxies, more 
than in any other bandpass. Our HDF-N/S clustering mea- 
surements are in agreement with the measurement from the 
much larger area of the WHDF. 

Our non-evolving models have some imp ortant differ- 



types -s e produce the other types by normalising the Sbc 
track to tne observed rest-trame colours. As m our earner 



ences with those used in the earl ier works of Roche et al 
(1993|) and [Metcalfe et al. (1995| ). Firstly, our models in- 



pap ers (Jones et al. 1991 ; Metcalfe et al. 1991; Metcalfe et al 
1995 ; Iv [cCracken et al. 200C ), the normalisations of our lu- 



minosity functions are chosen to match the galaxy counts at 
B ^ 18 — 20 and we seek to explain the low number counts 
at bright magnitudes from a combination of photometric er- 
rors and anomalous galaxy clustering, rather than substan- 
tial and hence unphysical evolution at low redshift in the 



elude the effects of internal extinction by dust ( correspond 
ing to Ab ~ 0.3 mag, using the dust model of Pci (1992)) 
and reddenin g by the Lyman alpha forest (as modelled in 
Madau (199'^ )). Both of these effects may become significant 
at the very faintest magnitudes we reach, where Zmed > 2. 
Secon dly, our fc— corrections ar e computed from the mod- 
els of Bruzual & Chariot (1993 ) for both our ev olvin g and 



luminosity of galaxies. Our models also mclude the ettects 



non-evolving models, whereas Roche et al. (1995) and Met 



of the Lyman-Qf forest, and, for spiral types, dust extinction 
correspon ding to the Large Magellanic Cloud as described in 
Pei (1992). The model redshift distributions produced are in 
good agreement w ith the redshift distributions of the CFRS 
(Lilly et al. 1995a) and from the Keck Hawaii redshift sur- 



calfe, Fong, & Shanks (1995) used polynom ial fits to the 



spectral energy distributions of Pence (1976) for their non- 
evolving models. These flts extend only to z ~ 2 and are 
held constant at higher redshifts. Thirdly, the redshift dis- 
tributions in these earlier papers were artificially truncated 
3. The sum effect of these differences is that in Roche 



at z ■ 



which tne inclusion of the evolutionary corrections have on 
our computed correlation function scaling relation, we also 
calculate an non-evolving redshift distribution. This is pro- 
duced by applying k— corrections only to each galaxy type. 



4 MEASURED AMPLITUDES 

In this Section we will present a comparison between our 
measurements of A^i and those in the literature. We defer 
an analysis of the implications these measurements have for 
the growth of galaxy clustering, as well as a discussion of 
our evolutionary models, to Section ^ here we present com- 
parisons only with the non-evolving, e = 0, qo — 0.05 model. 

In panels a-d of Figure ^ we plot our fitted correlation 
amplitudes extrapolated to one degree (filled symbols, cir- 
cles for WHDF and squares for HDF) as a function of sam- 
ple limiting magnitude for BRIK bandpasses in comparison 



vey (|C(|)wic, Songaila,^ fe Mu 199ti[j. lb illustrate the ettect gt al. (1993[) and [Metcalfe et al. (1995[) the slope of the 4, 



magnitude limit scaling relation remains constant whilst our 
slope begins to decrease at B ~ 26. By this magnitude limit 
the difference between our predictions and these previous 
works is ~ 0.2 in log(At^). 

Our 7?— band correlations plotted in panel (b) of Fig- 
ure [^ reach R < 26, although the number of galaxies in 
this catalogue is much smaller (~ 300) than in B— and 
consequently our errors are larger. Our measured clustering 
amplitud e at _R < 25.5 agrees well with t he faintest data 
point of Brainerd, Small, & Mould (1994); unfortunately. 



our survey area is too small to permit us to check our clus- 
tering amplitudes with values from the literature measured 
at brighte r magnitudes such as the large, ~ 2 deg^ CCD 
survey of Roche fc Eales (199i: ). Our measured clustering 
amplitudes in 7?— in the WHDF are much lower than the 
predictions of the non-evolving, stable clustering model. Our 
HDF clustering measurements are in good agreement with 
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Figure 3. The logarithm of the amplitude of the angular correlation function (jj{9) at one degree (A^) in the WHDF (filled circles), 
HDF-N (filled squares) and HDF-S (filled pentagons) shown as a function of apparent magnitude for BRIK selected samples (panels 
a-d) . For I, correlations axe plotted as a function of sample median magnitude. Open symbols show points from the literature. The solid 
line shows the predictions a non-evolving model with e = and with ro = 4.3/i~^ Mpc and qo = 0.05. Error bars on our measurements 
are calculated by a bootstrap resampling technique, as described in Section 3.1 . 
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the H DF clusterin g measurements of Villumsen, Freudling 
& Da Cfosta (199';^). 



cantly blu er by JC' ~ 22. The K— s elected (I—K) histograms 



For our /—band measurements, shown in panel (c) of 
Figure ^, we follow the practice in the literature and show 
correlation amplitudes as a function of sample median mag- 
nitudes and not limiting magnitudes. We follow the same 
proced ure for our model correlation amplitudes which are 



plotted at the median magnitude of each magnitude lim- 
ited sli ce. In addition to our J < 26 WHT data, we have a 



second, larger image taken at the INT which overlaps the 



WHDF. This covers a total of ~ 80 arcmin^ to / < 23.5 and 
allows us to determine A^i from Imed = 20 to Imed = 22 (the 
three brightest WHDF bins on the graph). The faintest bin 
in this INT dataset is in agreement with our measurements 
from the brightest bin of the WHT dataset. Furthermore, 
the preliminary result from the large-area 0.2 deg^ survey 
of Woods et al. (in preparation), shown as an open square, 
is agreement with our WHT measurement. At Imed ~ 26, 
measurement from the HDF fields appear to favour the 
lower values found in the WHDF. We note also that fainter 
Imed ~ 21, our measurement are below the predictions of 
the non-evolving e = model. 

Faintwards of Imed ~ 23 a discrepancy emerges be- 
tween our measurements and two previously published stud- 
ies. At Imed ~ 24, our WHDF clustering measureme nts are 



5 ti mes lower th an the measurements made by Brain 



erd & Small (1998) over two slightly smaller fields of area 
~ 30 arcmin at a similar limiting magnitude. At brighter 
magnitudes, our poi nts are also below the faintest bins of 
Postman et al. (1998). This work is a large-area CCD survey 



covering a contiguous 16 deg^ area and is currently the most 
reliable determination of galaxy clustering over wide angles 
and at intermediate {z ~ 1) depths. W e de fer a detailed 
analysis of these differences until Section 4.1 where we will 



attempt to quan tify if the discre pan cies between our s urvey 



and the works of Postman et al 



and 



Brainerd & Small could 



be explained in terms of cosmic variance effects. 

Finally, we turn to an investigation of galaxy correla- 
tions for K— selected samples. Until very recently measuring 
u{9) at near-infrared wavelengths was time-consuming and 
difficult as typical detectors covered only ~ 1 arcmin'^ . How- 
ever, wide-format IR arrays are becoming available making 
it now possible to conduct wide, deep surveys of the near- 
infrared sky. The filled circles in panel (d) of Figure ^ shows 
clustering amplitudes determined from our faint, H < 22.5, 
wide area (~ 50 arcmin^) Calar Alto Survey are shown, 
which will be described fully in a forthcoming paper (Mc- 
Cracken et al 2000, in preparation). Similarly, also plotted 
are clustering measurements from our 6' x 6' UKIRT IR- 
CAM3 mosaic ( [McCracken et al. 200o| ). At if ~ 27 we have 
computed a single point from NICMOS data taken as part of 
the north and south Hubble deep fields program (we trans- 
form from H to K using a model {H — K) colour). We note 
that all our measurements are in agreement with the predic- 
tions of our stable clustering, no luminosity evolution model. 

In plotting the H— limited Calar Alto points on our K— 
limited scaling relation we make two assumptions: firstly, at 
K ~ 22, {H — K) ~ 0.3; and secondly, for a given surface 
density, the clustering properties of H— selected and K— se- 
lected galaxies is identical. The first assumption seems rea- 
sonable, given that &t K ^ 20, galaxies in our survey have 
[H — K) ~ 0.3 and it is unlikely that they become signifi- 



shown in VlcCracken et al. (200C) support this. Also given 
that our Calar Alto H < 2Q agrees with our UKIRT 
K < 19.5 point, we conclude that our second assumption is 
also valid. 

Our points ai K = 19 — 20 agree with the survey 



Roche & Bales 



of Roche, Bales, & Hippelein (1998) and 
(1999D ; however at fainter magnitudes there is a dis crepancy 
between our amplitudes and the measurement of Carlber j 
et al. (1997). Once again, we defer a detailed discussion of 



the possible explanation of these differences until the follow- 
ing section. 

4.1 Quantifying errors in the correlation function 

In this Section we will investigate if we have estimated the 
magnitude of our correlation function error bars correctly. 
The small size of our field means our integral constraint 
(equation ^ corrections are large, and consequently accu- 
rate measurements of u!{0) axe dependent on an accurate de- 
termination of this quantity. Our main motivation is to see if 
we can explain the discrepancies between our m easurements 
of Au, at / < 25 and K ~ 21.5 with those of Brainerd &i 



Small (1998) and Carlberg et al. (1997). There are already 
indications that such "extra " variance could be si gnificant 
at the depths of our survey. Postman et al. (1998) directly 



address this question at shallower depths in their work which 
covers ~ 16 deg'^. By extracting 250 independent 16' x 16' 
fields from their survey (each of which is five times larger 
than the WHDF but at a brighter limiting magnitude) they 
find that the variance on (^(1') is comparable to its mean 
value of ~ 0.045, with extreme values reaching x3 this. Fur- 
thermore, they suggest that as the error distribution for A^, 
is non-Gaussian, and skewed positively, there could be many 
more areas in which A^j is below the mean value, rather than 
above it. 

To quantify the amount of "extra" variance which could 
affect clustering measurements in a very deep field like the 
WHDF we adopt a simple approa ch and generate large moc k 
catalogues using the method of Soneira & Peebles (1978). 



This is an purely empirical approach to generate a hierar- 
chically clustered distribution of points. We start by placing 
within a sphere of radius R a random distribution of sub- 
spheres each of radius R/\. Within each of these a further 
n spheres of radius RX^ are added. This continues through 
L levels; in our simulation we adopt L = 9. The amplitude 
of the correlation function is fixed by the number of centres 
used and the fraction of the total number of points which 
are retained; these quantities must be determined by trial 
and error. 

We measure the variance on the correlation function 
for many subsamples of this catalogue. We start by gen- 
erating a catalogue covering 6.25 deg^ with the same sur- 
face density of objects as in our real catalogue at / < 25 
(corresponding to ~ 7.5 x 10^ galaxies). Next, we mea- 
sure tjj{9) over the full simulated catalogue area. Our aim 
is to produce a catalogue for which the fitted correlation 
amplit ude Ior(Auj) at I < 25 is midway between the re- 
sult of 



Brainerd 



2-93Io.o6) and our 

own (log(Aa,) = -3.61^ „ 
sired correlation amplitude is produced it is randomly sub- 
sampled to produce 200 sub-areas each of which has the 
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Figure 4. Results from simulations of a 6.25 deg^ area with the same surface density of objects as our / < 25 catalogue. Panel (a) 
shows ijj(d) measured from an average of two hundred sub-areas each covering 50 arcmin^ (corresponding to the size of the WHDF), 
with error bars calculated using the normal bootstrap-resampling technique; panel (b) shows io{9) determined from the full simulation. 
Panels (c) and (d) show ^.'(6') as measured at 7 < 25 from the WHDF (again with error bars calculated using bootstrap-resampling) and 
the histogram of fitted values for \og{Ai^) from the simulations. The dotted lines shown on the histogram represent ztlcr deviations from 
the median value; the full simulation value is shown as the solid line and the average value as the dashed line. From this histogram we 
determine \og{A^) = -3.35l° j?- 



same field of view (~ 50 arcmin^) and galaxy surface den- 
sity at J ~ 25 as the WHDF (this translates to ~ 2000 
objects per field). On each of these sub-fields correlation 
amplitudes are measured using the same parameters as for 
the real data set, and a histogram is computed using each 
of these individual measurements of the simulated data. 

Figure ^ presents results from one set of simulations. 
Panel (b) in this Figure shows the measured correlation 
function for a synthetic catalogue generated using the 
method outlined above (error bars have not been plotted 
as they are smaller than the symbols for all bins). For this 
catalogue we find log(y4tj) = —3.26 for an integral con- 
straint C = 1. Panel (a) shows the average value of lo[6) 

from t^pn biTTirlrprl giih-gamplpg nf thra natalngiip ag ^vpll 



as the [fitted value to this average which we find to be 
log(A„) 



-3.32 



+0.08 
-0.09' 



For comparison, the fit to our / < 25 



observations is shown in panel (c); for the real data we find 
log(A.) = -3.6li°:i«. 



This procedure tests several important aspects of our 
technique. Because our simulated field is so large, the inte- 
gral constraint correction (Equation ^) which must be ap- 
plied to it is much smaller than the amount required for 
each of the individual sub-fields. Given that the A^^ which 
we measure from the full survey agrees to within the fit- 
ting errors to the A^^ determined from the average of two- 
hundred sub-fields we conclude that errors arising from an 
incorrect determination of the integral constraint are not 
significant. However, the agreement between the sub-fields 
and full-survey values is perhaps not surprising as both the 
catalogue and C were generated and calculated assuming a 
power-law slope 5 for uj{9) of —0.8. There is s ome indication 



that 5 becomes flatter at fainter magnitudes (Postman et al 



1996) although we are unable to test this with our current 
data set. A flatter correlation function at fainter magnitudes 
would lead to an underestimate of the integral constraint 
and a consequent underestimate of A^^. 
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More significantly for tliis current work, liowever, is tlie 
large dispersion we find for the fitted A^^'s from our sim- 
ulated fields. Panel (d) of Figure H illustrates this. From 
this diagram, log{A^) — — 3.35^q errors). This cor- 



responds to a linear error of ±2.2 x 10~ . By comparison, 
our errors d etermined from o ur bootstrap resamplings are 
±l.lx 10"''; Brainerd & Smail's errors are ±1.5x lO""*. Their 
points are an average of two widely separated ~ 30 arcmin^ 
fields and at J < 25 contain approximately the same num- 
bers of galaxies as our catalogue. On the basis of these simu- 
lations we conclude both our errors and those of Brainerd & 
Smailj^nderestimate the true error. Additional simulations 
at higher amplitudes also have higher variances; a second 
simulation at log(ylij) = —3.1 has a la error of ±2.6± 10~*. 
Adopting errors of this size, we find that our correlati on 
measurement is consistent with that of Brainerd & Small at 
the 2.5a level. 

Turning to the K- selected correlation amplitudes plot- 
ted in Figure ^ we no te that at K < 21.5 over ~ 27 arcmin^ 



population. As shown in Metcalfe et al. (200(]| ), the inte- 
grated star-formation history implied by our models agrees 
with current estimates of the global star-formation history 
of the Universe. The agreement between our predicted red- 
s hift distributions and observations at bright er magnitudes 



(Metcalfe et al. 1996; McCracken et al. 200C) gives us con- 



fidence in using these models to investigate the growth of 
galaxy clustering. 

In panels a-d of Figure ^ we plot our measurement of 
Aui for BRIK bandpasses (filled symbols) in comparison 
with the literature (open symbols), in addition to the pre- 
dictions of our evolutionary models with e = —1.2,0,1.0. 
We also consider a model with zero spatial curvature an d an 
non-monotonic ro — z relation, as explained in Section 3.2 



Carlberg et al. (1997) measure log(A 

at 



-2.72 



-1-0.08 
-0.10- 



This 



is also different from our work: at K < 21.7 we measure 
log(Atj) = — 3.53j;Q'36 in an area of ~ 44 arcmin'^. The num- 
ber density of galaxies at if ~ 21 is a pproximately th e same 
as at 7 ~ 25 



' 44 arcmin 

m 

and the amplitude of | Carlberg et al.'s point 



In comparing these four graphs, it is interesting to notice 
how the shape of the scaling relation is qualitatively differ- 
ent from bandpass to bandpass. To understand the origin 
of these differences, we start by emphasising that amplitude 
of the correlation function is directly related to the sam- 
ple median redshift. From Eq. (^) we see that A^ depends 
on the width of the redshift distribution, as well as its me- 
dian value. Our model scaling relations are therefore indica- 
tive of the median redshift of the model population. Other 
workers have commented on this relation previously ( Koo & 



is also within ~ 40% of Brainerd fc Smail's point. More- Szalay 1984 ); in this section we will attempt to see if our 



over, as we have described above, we find a larger variance 
on Lij{6) for simulations of higher amplitude. These consid- 
erations leads us to conclude that the stated error bar on 



Carlberg et al's measurement is also an underestimate of 
the true error. Furthermore, we conclude that our "low" re- 
sults at 7 ~ 25 and 7f ~ 21 are not inconsistent with the 
other results in the literature, given the large error bars af- 
fiicting measurements of A^^ in fields of this size and at these 
depths. 



5 INTERPRETATION AND DISCUSSION 
5.1 Galaxy clustering models 

After a decade of study, the general characteristics of galaxy 
evolution in the range z = Q — 1 have been broadly outlined, 
although many specific details have yet to be worked out, 
such as parameter dependence on morphology and intrinsic 
luminosity. Galaxy samples selected in bluer bandpasses are 
dominated by starburst populati ons as has been confirmed 
by many spectro.scopic .surypys I^LiUy, Cowie, OarH 



two-population luminosity evolution model can be used to 
explain the differences between the observed scaling of the 
correlation function amplitudes. 

Starting with the 7?— band scaling relation, we note that 
in this bandpass in all magnitude slices, the model popula- 
tion is dominated by spiral galaxies. Brighter than B ~ 22, 
evolutionary effects are negligible; however, faintwards of 
this, galaxies undergo 1 — 2 magnitudes of brightening, and a 
significant high-redshift tail becomes evident. For this rea- 
son the slope of the 73-band number counts is steepest in 
this range. The effect of this evolutionary brightening is 
to cause the median redshift of the B— selected redshift 
distributions to increase rapidly faintwards of 73 ~ 24; at 
B ~ 22, Zmed ~ 0.3, but by 73 ~ 24, Zmed ~ 0.7 and by 
B ~ 25, Zmed ~ 1- This causes A^^ to drop rapidly be- 
low the non-evolving prediction. This extended 73-band red- 
shift distribution, confirmed i n the spectroscopic survey of 



Cowie, Songaila, fc Hu (1996 ), allows us to explain the ob- 



served clustering amplitudes without recourse to positing a 
hypothetical, weakly clustered population domii iating the 



luer 



1991; Qlazebrook et al. 199£l Cowie. Songaila. fc Hu 1996 L 



B- selected samples, as did some earlier authors (Brainerd 



Small, fc Mould 1994; Infante fc Pritchet 1995a). In redder 



whereas samples selected in redder bandpasses show median 
redshifts and number count distributions which are closer to 



the non-evolving pred ictions (Vletcalfe et al. 1996 



Cowie, 



Songaila, fc Hu 1996| ). This mirrors broad trends seen in 



studies of the evolution of luminosity functions of colour- 



selected galaxy s amples (Lilly et al. 1995b ; Ellis et al. 1996 
Lin et al. 1997; Lin et al. 1999^, although distinguishing 



between density and luminosity evolution in these surveys 
is not straightforward. In interpreting our observations in 
terms of luminosity evolution models we are well aware of 
the limitations of these class of models, such as the difficulty 
in correctly reproducing the properties of near-infrared se- 
lected samples. Instead, we present our models as a sim- 
ple parametrisation of the observations using a model moti- 
vated by the star-formation behaviour of the entire galaxy 



bandpasses, the situation is slightly different; for example, 
in 7? spiral evolution is more gradual than in 73-, causing 
a much less pronounced slope change in the yli^-magnitude 
relation at 7? ~ 24. Similar considerations apply to the 7- 
band. By the 7f-band, however, the form of the Aui limiting 
magnitude relation is determined primarily by the early- 
type population; although the early-type counts turn over 
at ~ 20 they still comprise more than half of the to- 
tal galaxy population faintwards of this. Consequently, the 
A„-limiting magnitude relation has constant slope as these 
galaxy samples are dominated by slowly-evolving early-type 
populations. 
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Figure 5. The logarithm of the amplitude of the angular correlation function w(0) at one degree {Aa,) in the WHDF (filled circles), 
HDF-N (filled squares) and HDF-S (filled pentagons) shown as a function of apparent magnitude for BRIK selected samples (panels 
a-d) . For I, correlations are plotted as a function of sample median magnitude; open symbols show points from the literature. Error bars 
on our measurements are calculated by a bootstrap resampling technique as described in Section 3.1 . Also shown are the predictions of 
our best-fitting evolutionary model for three values of the clustering growth parameter e and for ro = 4.3/i~^ Mpc and go = 0.05. The 
long dashed line shows the predictions of the e = 0, qo = 0.5 dwarf-dominated model, and the dot-dashed line shows the predictions for 
the nA=0.7 case with dynamical evolution for haloes with rotation velocities >120 kms~^. 
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5.2 INIodel comparison with observations 



(20001 ), illustrated how the variations in IMF slope could 



From a visual inspection of Figure ^ our low-go evolutionary 
model with e = — 1 provides the best fit to the data at all 
magnitude limits; in the following Section we will present a 
quantitative analysis of the growth of clustering implied by 
our models. 

At fainter magnitudes, certain models are disfavoured; 
for example, the qo = 0.5, dwarf-dominated model, shown 
in Figure ^ as the long dashed line, produces much higher 
correlation amplitudes than the observations fainter than 
B ~ 26. This is because, in general, the median redshift of 
a magnitude-limited sample is lower for a low-go cosmology 
than for a high qo one, because the differential volume el- 
ement is smaller in the latter case. Consequently, for the 
same value of e and ro, A^i is higher for go ~ 0.5 than 
it is for go = 0.05. For this reason, our e = 0, go = 0.5 
dwarf model predicts higher clustering amplitudes than our 
standard e = go = 0.05 model. The exact magnitude of 
the differences between the two models depends (in addi- 
tion to the cosmological considerations outlined above) on 
where the median redshift of galaxies in the high-go model 
is greater than unity, where the star-formation rate for the 
dwarf types is constant, or less than unity, and where they 
rapidly fade. This high clustering amplitude leads us to re- 
ject the dwarf-dominated, e = 0, go = 0.5 model. Of course, 
this conclusion is dependent on the dwarf population having 
the same intrinsic clustering properties as the normal galaxy 



population, which may n ot be the case (Floche et al. 1993 
Infante fc Pritchet 1995a| ). 

Earlier works conducted in B— selected surveys sug- 
gest that by B ~ 26, A^i ceases to decline and reaches a 



cons tant, limiting value ( Roche et al. 1996 ; Metcalfe et al 
1995|; ^oche et al. 1993| ) ; more recently, [Brainerd fc Smai| 



(1998£]claimed to have detected a similar phenomenon at 
Imed ~ 24. The scaling relation for our _B— band evolution- 
ary model, shown in Figure |^ flattens off at very high num- 
ber densities (log(A'^9ai) ~ 6deg~^) and faint magnitudes 
{B ~ 28). At these limits, the relationship between number 
density and median redshift levels off. This is a consequence 
of the steep faint-end slope luminosity fun ction assumed fo r 
Scd and Sdm spiral galaxies (which have schechter (197C) 
function parameter a — —1.5), which means that at fainter 
magnitudes one observes intrinsically fainter rather than 
more distant galaxies, and also of the reduction of the cosmo- 
logical volume element at high redshift. Our _B— band corre- 
lations reach depths at which the correlation function is ex- 
pected to behave in this manner, and indeed from B = 27.0 
to i? = 28.0 we do observe that the amplitude of A^i is 
almost independent of magnitude. 

What are the implications of our K— selected A^i mea- 
surements? In McCracken et al. (2000) we demonstrated how 
the low me dian redshift found for K— sel ected redshift dis- 
tributions (Sowie, Songaila, fc Hu 199(^ ) placed stringent 



limits on the amount of evolution allowab le in these band- 
passes. In order for our PLE models to fit |Cowie, Songaila 



& Hujs] K < 19 redshift distribution (which has a very low 
median redshift, close to the predictions of a non-evolving 
model), we had to assume a steep slope {x = 3) for the initial 
mass function. This reduces the amount of passive evolution 
at K- for early-types, resulting in a total galaxv population 
with a lower median redshift (Figure 5 of McCracken et al 



affect the redshift distributions). 

As we can see from Figure ^ our e = 0, low go evolu- 
tionary model incorporating this steep IMF slope fits the ob- 
served clustering amplitudes for K— selected samples quite 
well. Therefore, the observed clustering amplitudes are con- 
sistent with the underlying redshift distribution for K- se- 
lected samples which has a low median redshift, close to the 
predictions of the non-evolving model. More significantly — 
and beyond the spectroscopic limit of the even the Keck 
telescope — our Calar Alto data at If ~ 22 indicates that 
even at these very faint magnitude levels, the H- selected 
galaxy correlations are still consistent with the non-evolving 
prediction. Galaxy merging, however, could provide another 
explanation for the low median redshift we infer for our K- 
sample. A low median redshift for K- selected surveys is a 
general prediction of the models of hierarchical galaxy for- 
mation (Kauffmarm fc Chariot 1998). 



5.3 Measuring the rate of clustering evolution 

In this Section we investigate what implications our mea- 
surements Aij have for the growth of galaxy clustering. As 
we have commented earlier, the small angular size of the 
WHDF means we are probing very small scales where the 
growth of galaxy clustering is expected to be highly non- 
linear. Most of the power in our correlation function signal 
comes from our inner bins, at angular scales of ~ 0.2'; at 
z ~ 1, the typical median redshift of our samples, this trans- 
lates to linear dimensions of ~ 0.05/i~^ Mpc (for go = 0.05). 
Additionally, how our samples are selected will affect clus- 
tering amplitudes. In our fiux-limited catalogues, a range of 
galaxy luminosities will be present, and local redshift sur- 
veys have shown that clustering amplitude may be a f unction 
of luminosi ty and morphology (Loveday et al. 1995; Tucker 



et al. 1997) 



With these caveats in mind, in Table 5.1 we present the 
results for best-fitting values for the parameters ro and e 
in Equation (^ determined by minimisation using our 
WHDF observa tion s and the model (A = 0, go = 0.05) out- 
lined in Section 3.2, As before, we use redshift distributions 
determined from our best-fitting evolutionary model. Be- 
cause of the strong co-variance between ro and e it is not 
possible to derive both parameters simultaneously from our 
dataset; instead we investigate what values of ro and e are 
implied by "reasonable" choices of these parameters. 

We wish to investigate what value of e best fits our 

data and to do this we fix ro to 4.3fe~^ Mpc. This value of 

ro is chosen to agree with angular correlation measurements 

determined from large Schmidt plate surveys. More recent 

work from local redshift surveys approximat ely agrees with 

this value. For example, Loveday et al. (1995) find for the bj 

selected APM an ro of 5.1±0.2/i"^ Mpc and 7 = 1.71. The 

R— selected Las Campanas Redshift Survey (Tucker et al 

1 



1993) finds ro = 5.0 ± 0.14/1" 
In general, we find e ^ 



Mpc. 
for ro 



4.3/1"^ Mpc and 
go ~ 0.05, from our own data alone. As we have already 
di scussed, in the I- band our points are different from those 
of Brainerd fc Small at the ~ 3a level. Their survey subtends 



^ 30' on the sky, and reaches similar depths to our own 
work, and so we would expect this survey to sample the same 
environments as our own, and therefore to show broadly 
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Bandpass e ro 

(ro = 4.3 Mpc) (e = 0) 



B 
R 
I 
K 



0.40 



+0.35 



3.70;^ 



2.65 



+1.30 
0.65 



2.30 



+0.75 
1.05 



1 lo+o'^s 

^•^^-0.50 



n QK+0.45 
•^•■^^-0.45 



0.05 



+0.65 

0.45 



4.30 



+0.70 
0.80 



Table 2. Best-fit values for e for ro = 4.3/i~^ Mpc and for rg for 
e = 0, using redshift distributions computed from our best-fitting 
evolutionary model and assuming go = 0.05. Errors quoted are 
ilcr. 
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Figure 6. Histogram of fitted values of log(Acj) carried out in 
100 sub-samples of a 1 deg^ area containing ~ 4 X 10® particles, 
representing the surface density of objects in our B ~ 27.5 sam- 
ple. The solid line shows the median value of the histogram, dot- 
dashed lines illustrate the ±lcr confidence limits, and the dashed 
line shows the fitted value obtained from the WHDF. 



similar gro-wth of clustering. Combining Brainerd & Small 's 
three I- limited points with our own, we deriwT^OTOT^, 
again for ro = 4.3/i~ Mpc. 

We have also ca rried out a simulation similar to those 
described in Section 4.1 to see how secure is our rejection 
of the co-moving amplitude in the range B — 27 — 28, the 
results of which are presented in Figure ^. This simulation 
has the same galaxy number density as our observations at 
B ~ 27.5 and covers an area of 1 deg^. It contains a total of 
~ 4 X 10® galaxies for an integral constraint C = 3.2. The 
simulation has \og{Aui) = —3.12 ± 0.02 (bootstrap errors), 
corresponding to the amplitude of our co-moving evolution- 
ary model at this magnitude. Errors calculated by resam- 



pling 100 WHDF-sized fields over this area gives a median 
log(A^) = -3.15tl:lt (la) and log(^^) = -3.151°;^^ (2a). 
At B ~ 27 we measure log{A^) — —3.7^0 j^s (bootstrap er- 
rors). Out of the 100 simulated fields, there are only two 
measurements at or below this value, leading us to conclude 
that in this magnitude range our measurement and the co- 
moving amplitude differs by at least 2a. 



5.4 The growth of clustering and biased galaxy 
formation 

How does this observed rate of clustering growth compare 
with measurements from the literature? With deeper pencil- 
beam spectroscopic surveys, it has become possible to mea- 
sure ro at successively earlier epochs and to use this to infer 
a value for e. Using a statistically complete subsamp le of 
591 galaxies fr om the Canada-Franc e Redshift Survey (Lilly 
et al. 1995a), LeFevre et al. (1996) were able to measure 
the evolution of ro in the interval < z < 1. They found 
ro(z = 0.53) = 1.5 ± 0.09 fe"^ Mpc (for go = 0.05), im- 



plying < e < 2. 
pie of 248 galaxies. 



Carlberg et al. (1997) 



found that for Mk < 



using a sam- 
-23.5 galaxies. 



ro(z ~ 0.6) = 2.0to'i Mpc. At higher redshift, [Carlberg 



et al. derive ro(z ~ 0.97) = l-4_o'2 h~ Mpc, which, com- 
bined with the lower redshift points from their survey, leads 
to e ~ 0.2 ± 0.5. The Canadian Network for Observational 
Cosmology (CNOC) have recently completed a large field 
galaxy survey i n the range < z < .7 with a sample size of 



10 galaxies (Carlberg et al. 1998). For luminous objects 



with corrected i?-band absolute magnitudes of M^'^ < — 20 
they find a slower clustering growth: e = —0.6 ± 0.4, with 
ro = 5.15 ± 0.15, strongly excluding clustering growth as 
rapid as e ~ 1. 

The large size of the errors on our fit does not permit 
us to make a detailed investigation of the dependence of e 
with sample selection. However, it is interesting that we find 
the rate of clustering growth to be slowest for our K- selected 
survey, and that our results are b roadly consistent with those 
from the K- selected sample of Carlberg et al. (1997). We 
expect our K- selected galaxies to be g ood tracers of the 
underlying matt er. A''- body simulations ( iGolin, Carlberg, &j 



Couchman 1997) have shown that in a low-density universe, 
the clustering of matter is expected to evolve as e ~ 0. 

Our finding that > e > 2, is in agreement with the 
expectations from biased models of galaxy formation, which 
find that at the < lh~^ Mpc scales we are sens itive to, clus- 
tering growth is relatively rapid (Figure 1. of Baugh et al 



(1999| ); Benson et al. (199^ )). In comparison, at larger scales 



(~ 5/i~ Mpc) the correlation function evolves much more 
slowly. At such separations, the clustering pattern is "frozen 
in" as the galaxies are tracing higher-mass haloes whose clus- 
tering evolution is close to e = —1.2. In Figure |^ we see that 
that the predictions from the non-zero lambda cosmology 
fitted to the growth of clustering as observed in the simula- 



tion of Colin, Klypin, 



Kravtsov (1999) is consistent with 
our observations and to the predictions of our e = model. 
The rapid decrease of the co-moving correlation length ro 
between z = Q and z ~ 2 for small haloes of V > 120km~^ 
is a prediction of biased models which find faster cluster- 
ing growth for intrinsically fainter galaxies. Even at brighter 
t han i? ~ 26 our sam ples are dominated by ~ L* galaxies 
( Metcalfe et al. 1996 ), unlike the Lyman-break galaxies of 
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Steidel et al. (1996) which higher intrinsic luminosities than 
this (The high resolution of Colin, Klypin, & Kravtsov 's sim- 



ulation and their adopted halo-finding algorithm makes it 
possible to locate haloes within haloes and therefore to suc- 
cessfully match the halo corr elation function with the APM 
galaxy correlation function (Maddox et al. 1990a). For this 



reeison we assume an approximate correspondence between 
these haloes and the intrinsically fainter galaxies which dom- 
inate our samples at B ~ 28.) 

Figure |^ indicates that our observations are consistent 
with models displaying the non-monotonic dip in the ro — z 
relation, and illustrates why the "epsilon" models have been 
so successful in describing the observed scaling relation of 
u!{0). If our models are correct, even to B ~ 28 the number 
of higher-redshift (2 > 2.5, more highly clustered galaxies 
forms only a small fraction of total sample size (less than 
<5%) and this explains why scaling relations in which the 
galaxy correlation length decreases monotonically z ~ 1 can 
successfully match the observations to B ~ 28. 



6 CONCLUSIONS AND SUMMARY 

In this paper we have presented a study of the projected 
two-point angular correlation function uj{8) as measured in 
the WHDF and compared our measurements to results ob- 
tained from the much smaller North and South Hubble Deep 
Fields. The clustering amplitudes determined from the HDF 
are consistent with those in the WHDF, but none of our con- 
clusions depend on our HDF measurements. In interpreting 
our results, we have used redshift distributions from a model 
which correctly predicts colours, number counts and redshift 
distributions for the faint galaxy population. 

We find that at a fixed separation the amplitude of w{9) 
measured in BRI bandpasses is lower than the predictions 
of a low-(/o model not containing luminosity evolution and 
in which clustering growth is stable in proper co-ordinates. 
For our K— selected samples, our correlation amplitudes 
are consistent with predictions from models having a low 
median redshift; we also find marginal evidence for a slower 
growth of clustering in these samples. 

If our evolutionary models provide a correct descrip- 
tion of the underlying redshift distributions (and compar- 
isons to available observations at brighter magnitudes sug- 
gest they do), then our WHDF clustering measurements 
are consistent with a clustering growth > e > 2 on the 
small scales (< lh~^ Mpc) which we probe. We have also 
shown that this result is consistent with prediction of bi- 
ased galaxy formation models which find faster clustering 
growth for intrinsically fainter galaxies like those which dom- 
inate our deep magnitude-limited surveys. We are able to use 
these rapid-growth "epsilon" models to successfully describe 
the clustering properties of our samples because the highly- 
clustered high-redshift galaxy population constitutes only a 
small fraction of the total galaxies observed in our survey. 

Finally, our constant correlation amplitude found at 
B ~ 27 — 28 is consistent with the expected reduction of 
cosmological volume element at high redshift and a steeper 
faint end slope for spiral galaxies, indicating that at these 
magnitude limits the median redshift of our sample ceases 
to increase. 
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